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Abstract 

We present an updated measurement of the time-dependent CP-violating asymmetry in B° — > 
K° S K° S K Q S decays based on 347 million T(4S) -> BB decays collected with the BABAR detector at 
the PEP-II asymmetric-energy B factory at SLAC. We obtain the CP asymmetries Sf = —0.66 ± 
0.26 ±0.08 and Cf = —0.14 ±0.22 ±0.05, where the first uncertainties are statistical and the second 
systematic. 
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1 Introduction 



In the Standard Model (SM) of particle physics, the decays B° — > K®K®K® are dominated by 
6 — > sss gluonic penguin amplitudes. Let 2/3 e // be the CP-violating phase difference between P° — > 
XgKgifg decays with and without mixing, and /3 = aig(—V c dV* h /VtdV* b ) where Vij are the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [1]. The difference | sin 2(3 — 
s'm2(3 e ff \ is expected to be nearly zero using calculations in SM, with theoretical uncertainties at 
the level of O(0. 01) [2], thanks to the factor \ V u bV us /VtbVt s \ which suppresses the other contributions 
in the SM. 

On the other hand, b — > sss decays involve one-loop transitions, so contributions from heavy 
new particles entering such loops can introduce new CP- violating phases, and these may contribute 
to P e ff [3]. The value of sin2/3 has been measured at the B factories in recent years with high 
precision [4, 5], with a world average of 0.685 ± 0.032 [6]. 

The Belle and BABAR collaborations have already reported measurements of CP asymmetries 
for B° -► (f>K s [7, 8, 9] (CP-odd) and B° K° S K° S K° S (CP-even) [10, 11]. 

The time-dependent CP asymmetry is obtained by measuring the proper-time difference At = 
tcp — *tag between a fully reconstructed decay B° — ► K®K®K® and the decay of a partially recon- 
structed tagging B meson (Ptag)- The expected asymmetry in the decay rate /+ (/_) when the 
tagging meson is a B° (B°) is given as 

e -\At\/r B() 

f±(At) = — x [1 ± S f sin (Am d At) =F Cf cos (Am d At) ] , 

4tbo 

where t b o is the B° lifetime and Amj is the B°-B° mixing frequency. The parameter Sf is non- 
zero if there is CP violation in the interference between decays with and without mixing, while a 
non-zero value for Cf would correspond to direct CP violation. In the limit of one dominant decay 
amplitude in b — > sss transition, the SM predicts no direct CP violation, and that Sf = — r/jsin2/3, 
within the theoretical hadronic uncertainties already mentioned. The CP eigenvalue r// = +1 for 
CP-even B° -► K S K° S K° S decays. 

In this paper, we update our measurement of the time-dependent CP-violating asymmetries in 
the decay B° — > K^K^R® previously presented in [11], using a larger data set and reconstructing 
the submode with one K® decaying into 7r°7r°. The absence of charged decay tracks originating at 
the B° decay vertex requires special techniques to deal with its reconstruction [12]. In addition, 
the final state has a definite CP content [13], so that an angular analysis is not needed. 

2 The BABAR detector and dataset 

The results presented here are based on 347.5 ± 3.8 million T{AS) — ► BB decays collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ collider, located at the Stanford Linear 
Accelerator Center. The BABAR detector [14] provides charged-particle tracking through a com- 
bination of a five-layer double-sided silicon microstrip detector (SVT) and a 40-layer central drift 
chamber, both operating in a 1.5 T magnetic field. Charged kaon and pion identification is achieved 
through measurements of particle energy-loss in the tracking system and Cherenkov cone angle in a 
detector of internally reflected Cherenkov light. A segmented CsI(Tl) electromagnetic calorimeter 
(EMC) provides photon detection and electron identification. Finally, the instrumented flux return 
of the magnet allows discrimination of muons from pions. 
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3 Bcp candidates selection 



The B° — > l£|if|.ffg candidate (Bqp) is reconstructed by combining three K® candidates. Two 
subsamples of Bqp candidates are formed. One subsample contains candidates formed by three 

Kg — > 7r + 7r - candidates in an event (B CP ^ )), while the other subsample is made by candidates 

formed by two K® — > ir + ir~ candidates and a third K® reconstructed in the 7r°7r° mode (-Bcp(oo))- 
The two subsamples have different signal to background ratios and therefore different analysis 
requirements were applied to obtain optimal selections. 

For the B CP ^ ) subsample we reconstruct K® — ► 7r + 7r~ candidates from pairs of oppositely 

charged tracks. The two-track composites must originate from a common vertex with a ir + ir~ 
invariant mass within 12 MeV/c 2 (about 4<r) of the nominal K® mass, and have a reconstructed 
flight distance (r^ ec ) between 0.2 and 40.0 cm from the beam spot in the plane transverse to the 
beam. We also require that the reconstructed K® has an angle between the transverse flight 
direction and the transverse momentum vector of less than 200mrad. For each Bqp/ + \ candi- 
date two nearly independent kinematic variables are computed, the beam-energy-substituted mass 

toes = \J ( s /2 + Pj • PB,) 2 /Ef + p B , and the energy difference AE = E* B - y/s/2. Here, (E^pj) is 
the four- vector of the initial e + e~ system, y/s is the center-of-mass energy, ps is the reconstructed 
momentum of the B° candidate, and E* B is its energy calculated in the e + e~ rest frame. For signal 
decays, the toes distribution peaks near the B° mass with a resolution of about 2.5 MeV/c 2 , and the 
AE distribution peaks near zero with a resolution of about 14MeV. We select B CP ( + _^ candidates 
within the window 5.22 < toes < 5.30 GeV/c 2 and -120 < AE < 120 MeV, which includes the 
signal peak and a "sideband" region for background characterization. 

For the i?cp(oo) subsample we form ir° — > 77 candidates from pairs of photon candidates in the 
EMC. Each photon is required to be isolated from any charged tracks, to carry a minimum energy 
of 50 MeV, and to have the expected lateral shower shape. We reconstruct K® — ► tt tt° candidates 
from 7T° pairs which form an invariant mass 480 < m n o n o < 520 MeV/ c 2 . -Bcp(oo) candidates are 
constrained to originate from the e + e~ interaction point using a geometric fit, based on a Kalman 
Filter [15]. We make a requirement on the consistency of the \ 2 of the fit which retains 93% of the 
signal events and rejects about 49% of other B decays. We extract the — ► ir + ir~ decay length 
L K o and the invariant mass (to 77 ) from this fit and require 100 < m 77 < 141 MeV/c 2 and L K o 
greater than 5 times its uncertainty. For K® — > ir + ir~ candidates we require 0.15 < r^ ec < 60.0 cm. 

For each B(jp(oo) candidate we compute two kinematic variables, the reconstructed mass m B 
and the missing mass m m i ss = \J (q e + e - — ?_b) 2 , where q e + e - is the four-momentum of the initial 
e + e _ system and q B is the four-momentum of the B° — > K®K®Kg candidate after a mass constraint 
on the is applied. By construction, the linear correlation coefficient between m m i ss and m B 
vanishes. This combination of variables shows smaller correlation (0.9% on reconstructed signal 
Monte Carlo events and 1.7% on the final data sample) and better background suppression with 

respect to the equivalent kinematic variables AE and toes used for Bcpu ) candidates. This 

is more relevant for -Bcp(oo) candidates given the asymmetric resolution on these variables due 
to 7T° energy reconstruction. We select -Bcp(oo) candidates within the window 5.11 < m m i ss < 
5.31 GeV/c 2 and —150 < m B —m B DG < 150 MeV/c 2 , where m B DG represents the nominal B° mass, 
reported by the Particle Data Group [16]. 

The sample of B° — > K®K®K® candidates is dominated by random K®K®K® combinations from 
e + e~ — > qq (q = u,d,s,c) fragmentation. Monte Carlo (MC) studies show that contributions from 
other B meson decays are small. We exploit topological observables to discriminate the jet-like 
e + e~ — > qq events from the more uniformly distributed BB events. In the T(AS) rest frame we 
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compute the angle 0% between the thrust axis of the B CP ^ ) (-Bcp(oo)) candidate and that of the 

remaining particles in the event. While |cos^| is highly peaked near 1 for e + e~ — > qq events, 
it is nearly uniformly distributed for BB events. We require |cos6*y| < 0.9 (0.95), reducing by 
one order of magnitude the number of background events. The maximum-likelihood fit described 
below also uses discriminant variables based on the momenta and angles of tracks in the event to 
discriminate Bqp candidates from qq. They are combined in a Fisher discriminant (J 7 ) [12] for 

Bq P u ) candidates, while in the case of -Bcp(oo) candidates we calculate the ratio L2/L0 of two 

Legendre monomials, defined as Lj = ^ i |p*|| cos#*p, where p* is the momentum of particle i in 
the e + e~ rest frame, 0* is the angle between p* and the thrust axis of the B candidate and the 
sum runs over all reconstructed particles except for the S-candidate daughters. 

After all selection requirements are applied, the average Bcp candidate multiplicity in events 
with at least one £>cp(q ) candidate is approximately 1.67, coming from multiple K® — > 7r°7r° 
combinations. In these events, we select the candidate with the smallest x 2 = — m K^) 2 1 a m i i 

where mi (m K o) is the measured (nominal K®) mass and a mi is the estimated uncertainty on the 
mass of the ith K® candidate. In simulated events, this selection criterion gives the right answer 
about 81% of the time. The remaining misreconstructed events, coming from fake K® — > 7r°7r° 
candidates, do not affect the determination of At and have a small impact on the other variables 

used in the final fit. The largest correlation is ~ 2.5%. In the case of events with B CP ^ ) 

candidates, only 1.4% of them have more than one candidate, and we apply the same criterion to 
select the best combination. 

Events coming from b — > ccs would reduce any sensitivity to departures from the Standard 
Model as this process is characterized by a Standard Model CP asymmetry (S ~ sin2/3 and C ~ 0). 

We therefore remove all B CP ^ ) (-Bcp(oo)) candidates that have a K®K® mass combination within 

3<7 (2a) of the XcO or Xc2 mass. After these vetoes the average efficiency, including K® sub branching 
fractions, is about 6% for B CP ^ + _y candidates and about 3% for B CP ( 00 y 

Combinatorics from other BB decays constitute a further source of background for -Bcp(oo) 
events. We take this into account by adding a component in the likelihood fit (see Sec. 5), where 
the shape of each likelihood variable is determined from a simulation of inclusive B decays. This 

contribution is found to be negligible in the case of B CP ^ ) events, and such a component is not 

included in the maximum likelihood fit. 

4 Flavor tagging and At reconstruction 

For each Bcp candidate we examine the remaining tracks in the event to determine the decay 
vertex position and the flavor of the -B^ag candidate. 

We use a neural network to determine the flavor of the B^ a g meson from kinematic and 
particle-identification information [17]. Each event is assigned to one of six mutually exclusive 
tagging categories, designed to combine flavor tags with similar performance and At resolution. 
We parameterize the performance of this algorithm with a data sample (Bf[ av ) of fully recon- 
structed B° — ► Z)M~7r + / p + ja\ decays. The effective tagging efficiency obtained from this sample 
isQ = J2 C e c (l — 2w c ) 2 = 0.305±0.004, where e c and w c are the efficiencies and mistag probabilities, 
respectively, for events tagged in category c. 

We compute the proper-time difference At = (zcp — zt ag )/7/3c using the known boost of the 
e + e~ system and the measured Az = zqp — ^tag, the difference of the reconstructed decay vertex 
positions of the Bqp and £>tag candidates along the boost direction (z). A description of the 
inclusive reconstruction of the &tag vertex is given in Ref. [18]. 
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To reconstruct the Bqp vertex from the K® trajectories we exploit the knowledge of the average 
interaction point (IP), which is determined from the spatial distribution of vertices from two-track 
events. We compute At and its uncertainty from a geometric fit to the T(4S) — > B°B° system 
that takes this IP constraint into account. We further improve the sensitivity to At by imposing 
a Gaussian constraint on the sum of the two B decay times (tcp + ^tag) to De ec L ua l to 2 t b o 
with an uncertainty of \[2 t b o, where t b o is the world average on the B° mean life [16], which 
effectively constrains the two vertices to be near the T(4S) line of flight [12]. The uncertainty on 
the IP position, which follows from the size of the interaction region, is on the order of 100 /im 
horizontally and roughly 4fim vertically. The mean uncertainty on zcp, a convolution of the 
interaction region and the vertex of the Bqp decay, is 75/um. The mean uncertainty on z tag is 
about 200/xm and thus the uncertainty in Az is dominated by the uncertainty in the vertex of the 
tagging decay. The resulting resolution is comparable to that in B° — > J/ip K® [12]. 

Simulation studies show that the procedure we use to determine the vertex for a Bqp decay 
provides an unbiased estimate of zcp- The estimate of the At error in an event reflects the strong 
dependence of the zqp resolution on the number of SVT layers traversed by the K® decay daughters. 
However, essentially all events have at least one K® candidate for which both tracks have at least 
one hit in the inner three SVT layers (at radii from 3.2 cm to 5.4 cm). In this case the mean At 
resolution is comparable to that in decays in which the vertex is directly reconstructed from charged 
particles originating at the B decay point [18]. For a small fraction (0.1%) of the signal events, 
at least one has tracks with hits in the outer two SVT layers (at radii 9.1cm to 14.4 cm) but 
none of the three K® s have hits in the inner three layers. In this case the resolution is nearly 
two times worse but the event can still be used in the CP fit. For events with a At > 2.5 ps or 
|At| > 20 ps, the At information is not used. However, since Cf can also be extracted from flavor 
tagging information alone, these events still contribute to the measurement of Cf. 

The At resolution function 1Z is parameterized as the sum of a 'core' and a 'tail' Gaussian 
distribution, each with a width and mean proportional to a At, and a third Gaussian with a mean 
of zero and a width fixed at 8ps [18]. We have verified on data that the parameters of 1Z for Bcp 
decays are similar to those obtained from the -Bfl av sample, even when the IP constrained vertexing 
technique is applied. Therefore, we extract these parameters from a fit to the -Bfl av sample. We find 
that the At distribution of background candidates is well described by a delta function convolved 
with a resolution function having the same functional form as that for the signal. The parameters 
of the background function are determined in the fit. 

5 Maximum Likelihood fit 

We extract the results from unbinned maximum-likelihood fits to the kinematic, event shape, and 
At variables. For each subsample we consider the logarithm of an extended likelihood function 

N T n 

c < N " v xU^.v/p;. 

* 3 

where Vj is the probability density function (PDF) for the ] th fit component and Nj the event 
yields of each component (N$ signal events, N q q qq events and, for -Bcp(oo) only, N BB BB decay 
events) Nt is the total number of events selected. The two C are then summed and maximized to 
determine the common Sf and Cf CP asymmetry parameters and the Nj which are specific to 
each subsample. 
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The At PDF for a given tagging category is V c (At, OAt)e c where e c is the tagging efficiency 
for tag category c. The total likelihood C is the product of likelihoods for each tagging category, 
and the free parameters are determined by maximizing the quantity ln£. Along with the CP 
asymmetries Sf and Cf, the fit extracts e c and other parameters for background. The background 
PDFs include parameters for the At-resolution function 1Z and for asymmetries in the rate of B° 
versus B° tags. We extract 43 parameters from the fit. 

The observables are sufficiently uncorrelated that we can construct the likelihoods as the prod- 
ucts of one-dimensional PDFs. The signal PDFs are parameterized from signal MC events. For 
background PDFs we determine the functional form from data in the sideband regions of the 
other observables where backgrounds dominate. We include these regions in the fitted sample and 
simultaneously extract the parameters of the background PDFs along with the fit results. 

There are 786 B CP ^ \ and 4550 B C p^ 00 ^ candidates that pass all the selection criteria. In 

Table 1 the events yields obtained in the fit are summarized for the two subsamples separately. 
Figure 1 shows the tties and AE (m m i ss and mg) distributions for B CP ^ + _^ (B CP ^) signal events 
with the sPlot event weighting technique [19]. The results of the fit are plotted as curves. 

6 Results 

The CP parameters extracted from the fit are summarized in Table 1. We obtain 

S f = -0.66 ±0.26 ±0.08, 
C f = -0.14 ±0.22 ±0.05, 

where the first error is statistical and the second systematic (evaluated as described in the following 
Section). The correlation between Sf and Cf is -8.5%. We evaluate the statistical significance 
of CP violation to be 2.Qa by calculating the 2AlogC variation when fitting data with Sf and 
Cf fixed to zero. Using a Monte Carlo technique, in which we assume that the measured values for 
the CP parameters on the combined data sample are the true values, we evaluate the probability of 
measuring the values reported in Table 1 for the two sub-samples. We find that the two sub-samples 
agree within 1.6er. 

Figure 2 shows distributions of At for i?°-tagged and S°-tagged events, and the asymmetry 
A(At) = {N B o — iVjjo) / (N B o + iVgo) > obtained with the sPlot event weighting technique [19]. 
This result is in good agreement with the value of sin2/3 from b — > ccs decays [6] . 





B CP(+-) 


Bcp(W) 


Combined 


N S 


116 ± 12 


60 ± 12 




N m 


670 ± 26 


4482 ± 71 




Nbb 




8 ± 25 




Sf 
Cf 


-1.04 +g;f 7 
-0.31 ±° ;i 


0.37 ± t u %i 
0.21 ± 0.38 


-0.66 ± 0.26 
-0.14 ± 0.22 



Table 1: Events yields and CP asymmetry parameters obtained in the fit. Statistical errors only 
are shown. 
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Figure 2: Distributions of At for events weighted with the sPlot technique for B tag tagged as B° 
(top) or B° (center), and the asymmetry A(At) (bottom). The points are weighted data and the 
curves are the corresponding PDF projections. 
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7 Systematic studies 



We obtain systematic uncertainties in the CP coefficients Sf and Cf due to the parameterization 
of kinematic variables and event shape PDFs in signal and background by varying the parameters 
within one standard deviation (evaluated from a fit to Monte Carlo simulated events) . There might 
be a contribution to Cf and Sf from CP violation in the BB background. In the fit, the values 
of the effective CP parameters (S BB and C BB ) for the BB background are fixed to zero. They are 
varied within the whole physical allowed range S BB + C 2 BB < 1 and we take the largest variation 
on signal Sf and Cf as systematic uncertainty. 

We evaluate the uncertainties associated with the assumed parameterization of the At resolution 
function for signal and .B-B-background by varying the parameters within one standard deviation 
(extracted from a fit to the -Bfl av sample). The uncertainties due to knowledge of efficiencies 
and dilutions of flavor tagging and possible difference in the efficiency between B° and B° are 
evaluated in the same way. The mass difference between the two B° mass eigenstates, Arrid, and 
the B° mean life, t b o , values held fixed in the fit are varied within their uncertainties determined 
in world averages [16]. 

We also estimate different uncertainties associated with vertexing. We take the largest value 
of Sf(Cf)st — S f(Cf )true from fits to signal Monte Carlo events where realistic misalignments of 
the SVT silicon wafers have been introduced. Here the Sf(Cf)&t represents the result of the fit to 
these simulated events, while Sf(Cf)t Tue represents input values in the Monte Carlo generation. We 
include an additional contribution from the comparison of the description of the resolution function 
(RF) between IP-constrained vertexing and nominal vertexing in the case of B° — ► J/ipKg events. 

We assign a systematic uncertainty on our knowledge of the beam spot position by shifting the 
beam position in the simulation by ±20 jim in the vertical direction. The sensitivity due to any 
calibration problems or time-dependent effects is evaluated by smearing the beam-spot position 
by an additional ±20 jum in the vertical direction. The effect of neglecting possible correlations 
between the variables in the fit is estimated with a Monte Carlo technique. We also estimate the 
errors due to the effect of doubly CKM-suppressed decays on the tag side by varying the value of 
the rate of such decays and the strong and weak phase within conservative limits [20] . 

We add these contributions in quadrature to obtain the total systematic uncertainty. The 
summary is reported in Table 2. The largest contributions are related to the knowledge of the PDF 
parameters and the CP content of the BB background. 

8 Conclusions 

In summary, we have measured the time-dependent CP-violating asymmetries for the decay B° — > 
K° S K° S K° S : Sf = -0.66 ± 0.26 ± 0.08 and C f = -0.14 ± 0.22 ± 0.05. Within the current experimen- 
tal uncertainties, these measurements are in good agreement with the SM expectation. 
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data TCF 
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0.003 
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Table 2: Summary of systematic uncertainties on S and C. 
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Abstract 



We present an updated measurement of the time-dependent CP-violating 
asymmetry in B° — » KgKgKg decays based on 347 million T(4S) — > -BP 
decays collected with the BABAR detector at the PEP-II asymmetric-energy 
B factory at SLAC. We obtain the CP asymmetries S f = -0.66±0.26±0.08 
and Cf = —0.14 ± 0.22 ± 0.05, where the first uncertainties are statistical 
and the second systematic. 
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